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Abstract 



Self-attracting walks {SATW) with attractive interaction n > display a 
swelling-collapse transition at a critical Uc for dimensions d>2, analogous to 
the Q transition of polymers. We are interested in the structure of the clusters 
generated by SATW below Uc (swollen walk), above Uc (collapsed walk), and 
at Uc, which can be characterized by the fractal dimensions of the clusters df 
and their interface di. Using scaling arguments and Monte Carlo simulations, 
we find that for u < Uc, the structures are in the universality class of clusters 
generated by simple random walks. For u > Uc, the clusters are compact, 
i.e. d{ = d and di = d — 1. At u^, the SATW is in a new universality class. 
The clusters are compact in both d = 2 and d = 3, but their interface is 
fractal: di = 1.50 ± 0.01 and 2.73 ± 0.03 in d = 2 and d = 3, respectively. In 
d = 1, where the walk is collapsed for all u and no swelling-collapse transition 
exists, we derive analytical expressions for the average number of visited sites 
(S) and the mean time {t) to visit S sites. 

PACS numbers: 68.35.Rh, 64.60.Fr, 05.40. -a 



1 



Typeset using REVTgX 



I. INTRODUCTION 



In previous years various models of random walks (RW) with memory or interaction 
have been studied p|-p!0|, in order to account for distinct features of physical, chemical, and 



biological systems whose complexity goes beyond what can be obtained from the simple 
random walk picture. Perhaps the most extensive studied model is the self-avoiding walk 
(SAW), where the random walker is not permitted to step on already visited sites, simulating 
the behavior of linear polymers. Although all investigated RW models with memory are 
similar in the sense that they incorporate interactions between steps, they display quite 
distinct asymptotic properties. Therefore, they belong to universality classes which are 
usually different from RW or SAW, though they may cross over to either RW or S'ylVK 
behavior in some limits. Common properties for describing the behavior of a walker are the 
exponents u and k, characterizing the scaling with time t of the mean square end-to-end 
distance 

{R\t)) ~ (la) 
and the average number of visited sites 

{S{t))r^t\ (lb) 

For RW the exponents are z/ = 1/2 for all dimensions d, k = 1/2 for d = 1 and k = 1 for 



d>2^. For SAW, A; = 1 for all d and = 3/{d + 2) [y,y]. A comparative study ITj of 
some of these models p|,|5|,p|JT0|] in one dimension has shown that the characteristic exponents 
depend crucially on the particular form of the interaction between the steps. Some of the 
important mechanisms are the range of the interaction, the presence of cumulative memory 
effects, and the global or local normalization conditions. Models with global and local 
normalization conditions are also refered to in the literature as static and dynamic models, 
respectively. 

Recently, one of the dynamic models without cumulative memory effects, the self- 
attracting random walk (SATW) |0, has been found to display, in contrast to all other 



models, a swelling-collapse transition at a critical attractive interaction in (j > 2 [[15 
The characteristic exponents u and k for SATW are in different universality classes for be- 
low Uc (swollen walk), at Uc, and above Uc (collapsed walk). Above and below criticality, u 
and k have been determined analytically. At criticality, the exponents could only be studied 
numerically, and due to the finite-size effects close to a transition, there have remained open 
questions regarding the asymptotic behavior of SATW. A careful analysis of the simulation 
data and a scaling approach different from the one developed in is found to be neces- 
sary for a comprehensive study of the structural properties of SATW, especially at critical 
u = Uc, which is the focus of this paper. To determine the fractal dimension of the cluster 
and the interface generated by the walker, and to give more precise results for the charac- 
teristic exponents, we investigate the temporal development of the number of sites visited 
by the walker which have a certain amount of already visited next nearest neighbor sites. 
Identifying the sites belonging to the external and internal perimeter of the cluster in ci = 2, 
the fractal dimensions of these structures are studied for attractions below, at, and above 
criticality. 

The paper is organized as follows: In Sec. |T| several known RW models with interactions 
are summarized and their behavior is compared. The SATW model is briefly reviewed 
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in Sec. |ITT[ presenting the analytical and numerical results for the mean square end-to-end 
distance and the average number of visited sites for varying strength of attractive interaction 
following [|T5| , showing the evidence for the phase transition. In Sec. |V| we investigate the 
structure of the cluster grown by the walker for different attractive interactions using a new 
approach consisting of scaling arguments and Monte Carlo simulations which leads to more 
insight for the behavior of the system at criticality. We study the fractal dimensions of the 
cluster and its interface in c? = 2 and d = 3. Closed form expressions for the average number 
of visited sites and the mean time to visit a certain number of sites in d = 1 are given in 
Sec. |V[ Finally, in Sec. |V| we summarize our main results. 



II. RW MODELS WITH INTERACTION 

When a self-avoiding constraint like in the S'ydVK model is added to the RW, the evolution 
of the walk becomes heavily dependent on the entire history of the walk, converting it into 
a non-Markovian system. A bridge between ordinary RW and non-Markovian walks can 
be constructed by associating energies Etiw) to all possible configurations w of a t step 
random walk, defining an ensemble probability for a certain walk configuration as Pt{w) = 
exp[—PEt{w)]/'^^exp[—(3Et{w)], with P = (k^T)^^, where /cb is Boltzmann's constant 
and T is the absolute temperature. In the high-temperature limit all walks become equally 
likely, but at finite temperatures the ensemble probabilities of individual paths differ. If 
Et{w) < for walks with many visited sites, walks prefering to explore new terrain dominate 
the system, whereas if Et{w) > for walks with many visited sites, the system is governed 
by configurations of walks attracted to their own path. Models based on this concept are 
by definition in the class of static models, for a comprehensive overview see One 
straightforward way of modelling attractive and repulsive interaction between steps, known 
as the interacting walk [|l|, is to assign an energy Et{w) = g St{w) to a walk configuration, 
where St{w) is the number of visited sites of a t step RW configuration. For interaction 
parameter g = the simple random walk is recovered, while for g < the walk becomes 
repulsive, and the characteristic exponents are i> = k = lind=l 0,^. For g > the 
walk is attractive and k = 1/3 and 2/3, in c? = 1 and 2, respectively, as well as d/ {d + 2) in 



d>3 mUM 



In a more generalized static model including cumulative memory effects the energy of 
a walk configuration is Et{w, a) = g "^f (s), with < a < 2, where nt{s) is the number of 
times a certain site s has been visited after t steps. For cumulative memory parameter a = 
the previous model of is recovered, and for a = 2 it is also known as the Domb- Joyce 
model while for (7 = or a = 1 it reduces to simple random walks. The model is 

repulsive for a < 1 and (? < or a > 1 and (7 > 0, in the latter case Flory arguments give 
V = {a + 1) / {2 + ad — d) independent of g for d < dc = 2a/(a — 1), where dc is the critical 
dimension. For the attractive regime results are only known in d = 1, where for a < 1 and 
(7 > it exhibits continuous varying exponents depending on a with u = k = (1 — a)/(3 — a), 
while for the other attractive branch a > 1 and (7 < it is always self-trapping as (R) and 
(S) saturate. 

An approach analogous to the one above for the static models can be made for the 
less investigated dynamic models, where the local normalization is done by assigning a 
probability Pj to the walker to step to the next nearest neighbor site i during the evolution 
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of the walk, with Pi = exp{g / Y^f^^exp{g nf), a > 0. Here, is the number of times 
the neighbor site i has aheady been visited in the previous t time steps For g = the 
simple random walk is recovered, while for g < the walk is repulsive and the exponent u 
is determined to be z/ = {a + l)/{2 + a d) for d < dc = 2. For the attractive regime > it 
is known that in c? = 1 the walk is always self-trapping. A special case of this model is the 
true self- avoiding walk of with g < and a = 1. 

In surprising contrast to the results for all above mentioned models, where the character- 
istic exponents are always independent of the actual strength of the attraction or repulsion 
parameter g, are the results for the SATW model focused on here. In this model a 
random walker jumps with probability P ~ exp{un) |16| to a next nearest neighbor site, 
with n = 1 for already visited sites and n = for not visited sites. For m > 0, the walk 
is attracted to its own trajectory, so that SATW is an extension of with attraction pa- 
rameter u = g > in the limit of no cumulative memory effect a — 0. Note that the 
results of Ref. can not be directly apphed to the SATW model as they are based on 
the restriction a > 0. Some representative examples of two dimensional clusters grown by 
SATW for different values of u at three distinct times t of evolution are shown in Fig. |l|. 
The exponents k and u of SATW have been found to depend on u |]^|T^-|T^|, although it 



was not clear for some time if u and k decrease continuously with increasing u |18,19|, or if 



a critical value Uc exists, below, at, and above which the exponents show different universal 
behavior • Recently, it was found by exhaustive computer simulations up to t = 5 ■ 10^ 
time steps that there exists a swelling- collapse transition for SATW at a critical attraction 
Uc W^, analogous to the B transition in linear polymers at temperature T = O when an 



attraction term exp(— A/T), A < 0, is added to the self-avoiding constraint [|T2|,|13|. This 



phenomenon of a swelling-collapse transition can only occur because of a balance of the 
interaction energy and the configurational entropy of the SATW at criticality. When the 
energy is not of the order of the entropy as investigated in for a > 0, the walk collapses 
for any attraction u = g > 0. 



III. SWELLING-COLLAPSE TRANSITION FOR SATW 



In [|T5| it was shown that the characteristic exponents u and k for SATW are in different 



universality classes for u < Uc, u = Uc, and u > Uc- Above Uc, the walk collapses for d > 2, 
and u and k are given by 

^ (2a) 



d + 1 
and 

Eqs. (H) follow (cf. fl^), since for sufficiently strong attraction u > Uc the grown clus- 
ters are compact (see Fig.|I]), so the average number of visited sites scales with the root 
mean square displacement (i?) = (i?^)^^^ as (5) ~ {R^ ■ Comparing this to Eqs. (|l]) 
yields k = vd for u > Uc- Also, the mean cluster growth rate is proportional to the 
ratio between the number of boundary sites and the total number of the cluster sites, 
d{S)/dt ~ {Rf-^/{Rf ~ t-" I^Jig]. Thus {S) ~ t"'+\ Combining these results with 
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Eq. dUb), one obtains Eqs. (0). Below a critical interaction Uc, the walk swells and the 



exponents are as with no interaction [11|, i.e. u = 1/2 and k = 1 for d > 2. The above 
analytic arguments are in good agreement with numerical simulations in d = 2 and d = 3 
(see Fig. |^). At Uc, the exponents are numerically determined to be = 0.40 ± 0.01 and 



kr = 0.80 ± 0.01 for d = 2 as well as z/, = 0.32 ± 0.01 and fc, = 0.91 ± 0.03 for ci = 3 [15 



Note that for = 1 no swelling-collapse transition exists, as the walk is collapsed for all u, 



and Eqs. reveals the known values k = 1/2 and u = 1/2 for random walk in c? = 1 |TT|pT 
(see also Sec. 0). 

In the analogous static model of the interacting walk [|I| with attractive interaction and 
no cumulative memory effects, a phase transition can not be observed due to the fact that 
the global normalization condition increases the weight of the interaction energy more than 



the local normalization condition, as already observed in for d = 1. This can also be seen 
from the asymptotic behavior of the characteristic exponent k = d/{d + 2) for the attractive 
interacting walk [|T],|TT| in comparison to the 'less collapsed' k = d/{d+ 1), Eq. (^), for 
the SATW discussed here. Therefore, a static interacting walk with the slightest attraction 
u = g > has a qualitatively different behavior than the ordinary RW and is collapsed. Note 
that in both models k never becomes independent of the dimension, although it approaches 
unity from below in the large dimensionality limit. 

Due to finite time effects at criticality it is not possible to determine the exponents Uc 
and kc more and more accurately by simply increasing the number t of time steps performed 
by the walker. As long as the attraction is slightly above or below criticality, asit is always 
the case in numerical simulations, the exponents will finally cross over to their expected 
values above resp. below Uc after some time t. Introducing a crossover time t^, below which 
the exponent u is close to and above which u approaches 1/2 for u < and l/{d+l) for 



u > Uc (see Fig. 1 of |]TS[), the following scaling theory holds: 

(i?(t))~r^/±(t/t^), (3a) 

with 

= \u — Uc\~°' , (3b) 



where the plus sign refers to u > Uc, the minus sign to u < Uc, and the exponent a has 
to be determined numerically. As is the only relevant time scale, the scaling functions 
bridge the short time and the long time regime. To match both regimes, it is require that 
f±{x) = const for X < 1 (t -C t^), and /+(x) ~ a;i/(<i+i)-!^c^ ^ a;i/2-!^c fg^. 2; > i. An 

analogous scaling approach holds for {S{t)), and an excellent data collapse can be obtained 
for a = 7±lmd = 2 and a = 5.0 ± 0.5 in ci = 3 (see Fig. 3 of [|T3[), confirming the 
numerical values for Uc, kc, and Uc determined from Fig. |^. 

Since the mass of the cluster generated by the walker, consisting of all visited sites, scales 

as 

is) - (4) 

the ratio k/u corresponds to the fractal dimension d{ of the cluster, 

d, = ^. (5) 
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In (i = 2 the clusters are compact for aA\ u as k / u = d{ = d. In d = 3 they are compact for 
u > Uc, while for u < Uc, the fractal dimension of clusters generated by simple random walks 
d{ = 2 < d is obtained. At the criticality, d{ was found numerically to be df = 2.84 ± 0.25, 
but the possibility that df = d could not be ruled out in Ref. ||T5[| . 



IV. THE CLUSTER AND ITS INTERFACE AT CRITICALITY 

To clarify if the grown clusters in d = 3 at u = Uc are compact or fractal, and to 
learn more about the structure of the .SylTVF-clusters and their interfaces at criticality, we 
consider the following: Denoting by Ni{t) the number of cluster sites which have i of their 
2 d next nearest neighbor sites belonging to the cluster after t time steps, then the number 
of all cluster sites S(t) is the sum of all Ni{t) 

S{t) = N,{t) + N2{t) + ... + N2d{t). (6) 

The cluster growth rate is equal to the probability to be on the boundary of the cluster 
multiplied by the conditional probability to expand the cluster while being on its boundary. 
Suppose the walker is on a site which has i visited next nearest neighbor sites. As the 
probability to step to a next nearest neighbor site is P ~ exp(Mn), with n = 1 for already 
visited sites and = for unvisited sites, the probability to jump on a visited neighbor site 
is proportional to i exp(M), whereas the probability to jump on an unvisited neighbor site is 
proportional to 2d — i. Thus the normalized probability Pj to expand, i.e. to jump on one 
of the 2d — i unvisited next nearest neighbor sites, is given by 



2d 



i exp(-u) + 2d 
Therefore, the average cluster growth rate is 



(7) 



The situation i = 2d is special because P2d = 0, as the cluster cannot be expanded from 
a site where all surrounding sites already belong to the cluster. Assuming that the average 
number of cluster sites (Ni) which have i already visited nearest neighbor sites scales as 

{Ni) ~ {Sf' with < < 1, (9) 

there must be at least one i for which a, = 1 to ensure that Eq. (Q) holds. In general, one 
can distinguish between two different cases: 

(i) there exists at least one i < 2d ioT which = 1 

(ii) only = 1. 

In case (i), the average cluster growth rate d{S)/dt is dominated by the (Ni) for which 
ai = 1, leading to d{S)/dt ~ Pi {N,)/{S) = Pi (S)"V(S) = const. Therefore, in this 
case one gets k = 1 as {S) ~ t. In case (ii), when only 02 d = 1, the cluster growth 



6 



rate is dominated by the (Ni) for which = Omax = max(ai, a2, a2d-i)- This gives 
d{S)/dt ~ (^)"™/(^) = and hence 

^ "max 

when compared to Eq. (|l]b). As shown below, the grown cluster is fractal for case (i) while 
for case (ii) it is compact. 

The above considerations also enable us to determine numerically the exponent A; in a 



way different to Ref. ||T^, and additionally to gain more insight into the grown structures, 
especially in the regime at u = Uc- Fig. || shows that the assumption of Eq. (|^) is clearly 
supported by numerical results. When plotting {Ni) versus {S), the decision to which case 
(i) or (ii) a certain regime belongs becomes obvious from Fig. ^ for all regimes except for 
u = Uc in d = 3. For -u < Mc in (i = 3 we observe case (i) (see Fig. |^), as = 1 for all 
i. The other regimes belong to case (ii) (see Figs, ^a, and ^), since only a2d = 1. 

At u = tic in c? = 3 (Fig. |^), a more detailed investigation of the results is necessary. By 
plotting the successive slopes Oj = dln(A^i)/d In(S') of the data of Fig. |e versus l/ln(S') 
(Fig. |a), it becomes obvious that only a2d = 1- Thus this regime is belonging to case (ii). 
Moreover, Fig. ^ indicates that all ai,a2, 02^-1 asymptotically have the same value, so 
c^max = Oi = 02 = ... = 02^-1. As this is also found for all other regimes, in the following we 
will denote oj = Oj for all 1 < i < 2 d — 1 in both cases (i) and (ii), in distinction to a2d- The 
values of the exponents are summarized in Table confirming the previous results for the 
exponent k (see Table |) when following Eq. ([T0|). At u = Uc iri d = 3, kc can be determined 
more precisely than in Ref. [|T^] to be k^ = 0.91 ± 0.01. Note that another approach to 
confirm our numerical results for k is to plot (A^,) versus t in a double logarithmic plot. 
Denoting the resulting slopes as Kj, we expect 

Ki = kai (11) 

as (Ni) ~ (S)^' t^"-^ ~ t'^^ when combining Eqs. (|l]b) and (P). The numerical values 
determined for Ki are summarized in Table y and are in excellent agreement with the values 
determined for when comparing with Eq. (pTj). 

When examining the distinction between cases (i) and (ii), it becomes clear that the 
structures grown by the walker in case (i) are fractal, while the ones in case (ii) are compact. 
This can be explained by considering that in case (ii) the sites which do not belong to the 
interface dominate the growth process as a2d > ai, leading to a compact structure (i.e. 
d{ = d). In case (i) all sites have the same contribution to the growth of the structure due 
to a2d = CLi [0]- As we found that u = Uc in d = 3 belongs to case (ii) (see Fig. ||a), 
the structure must be compact, and we can definitely conclude that in this regime d{ has 
to be (if = d = 3, correcting the value 2.84 ± 0.25 obtained earlier [0. The latter value 
was calculated by combining the numerical results for the exponents kc = 0.91 ± 0.03 and 
z/c = 0.32 ±0.01 following Eq. (|^). From the new numerical result k^ = 0.91 ±0.01 and 
^'c = kc/df = kc/d, we obtain the more accurate estimate = 0.303 ± 0.005. 

Next we focus on the interface of the grown cluster. The total interface can be divided 
into the external perimeter, which is usually the more interesting fraction as it constitutes 
the reaction front with the environment, and the internal perimeter, which is the boundary 
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of the inner holes of the cluster p4|. First we investigate the total interface. Its mass / is 
equal to Ni + N2 + ■■■ + A^2d-i and therefore scales as 



Combining Eqs. (^) and 
di of the interface is simply 



(/) ~ {sr. (12) 

with Eq. (|1^) yields (/) ~ {R)"'^'^^ . Hence the fractal dimension 



dj = ai d{ 



(13) 



The values for di according to Eq. ( ]T3| ) are given in Table fT|. For u < Uc, we recover the 
values for normal random walks p] , p5| , di = 2 in d = 2 and 3. As the fractal dimension 
of the external perimeter of a random walk is c/ep = 4/3 < rfi in d 
that the internal perimeter is dominating the interface for u < Uc- In contrast, in d 

m 



this suggests 
3 the 



external perimeter governs the interface, dEp = 2 = di 



as three dimensional holes are 
less likely than two dimensional ones due to geometrical constraints. For u > Uc, we find 
di = d — 1, clearly confirming the assumption that the structure of this regime is collapsed 
(see Fig. |l]). It forms a compact disk resp. a compact sphere with a rather smooth surface 
and holes only in the surface layer. The clusters are similar to the clusters grown in the 
Eden model where each unvisited next nearest neighbor site of the cluster has the same 
probability to be occupied at the given time step. It remains an open question whether 
these clusters are in the same universality class, as one has to check if the surface for u> 
resembles the self-afiine surface of Eden cluster. At u = Uc, the structure of the interface is 
fractal in both d = 2 and 3 with the novel values of di = 1.50 ± 0.01 and di = 2.73 ± 0.03, 
respectively. The rather surprising result for di in d = 3 can be understood when considering 
that the structure of the interface at m = Mc is more compact than for u < but the cluster 
is not collapsed as for u> u^. 

A straightforward way to confirm the above numerical results for di is to determine the 
fractal dimension of the interface by measuring it directly on the grown clusters instead of 
calculating it dynamically during the growth process. Using the sandbox method |2^, the 
outcome supports the results stated above, although they are less precise due to the fact 
that the information about the growth process is lost. For the results of di at criticality 
in d = 2 and d = 3 see Fig. |^, showing a good agreement with the values summarized in 
Table |T| and clearly excluding the value of d{ = 2.84 determined in for d = 3, which 



would lead to di = 2.58 when following Eq. (|T3D. 

A further interesting question is which perimeter, external or internal, is dominating the 
interface of the clusters at criticality, and what is the actual value of the non-dominating 
perimeter. Using the approach of Ref. to identify the sites A^ep belonging to the external 
perimeter of a cluster in d = 2, and identifying the internal perimeter sites Njp as the 
interface sites not belonging to A^ep, we expect the average number of external perimeter 
sites (A'ep) to scale as 



and the average number of internal perimeter sites (Ajp) to scale as 



(14a) 
(14b) 



Following the reasoning leading to Eq. (|T3D yields 
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c^EP = Qep di (15a) 

and 

dip = aip c?f. (15b) 

The results for oep (see Fig. ||) recover the expected behavior of d^p = 4/3 for m < Uc, and 
indicate that at criticahty d^p is clearly smaller than for below Uc, although it is difficult to 
determine the precise value of c^ep from this plot. Thus, at Uc the interface is dominated by 
the internal perimeter as c^ep < di. This is confirmed by the results for aip at criticality, were 
the expected value of dip = rfj = 1.5 is recovered, although with unsatisfactory precision p9| . 



Using the sandbox method [^] to measure the fractal dimension of the perimeter directly, 
the results support dip = 1.5 more precisely, and we find d^p = 1.25 ± 0.05 for u = Uc (see 
Fig. 1^). The fact that c^ep for u < Ucis definitely smaller than c^ep at m = Mc reinforces that 
at criticality the SATW model is in a new universality class. For u > Uc one might expect 
that the interface is dominated by the external perimeter due to the collapsed structure of 
the cluster. However, we find that the internal and the external perimeter contribute equally 
to the interface, dip = dEp = I = di (see Fig. ^). 

In d = 3, the fractal dimension of the external perimeter is c^ep = 2 = di for u < Uc, 



supposing that c^ep follows the random walk result [TT[] like all other characteristic values 
determined for the SATW in this regime. Unfortunately it is not possible to determine c^ep 
or dip in d = 3 numerically using the algorithm of ||2^ to check which one, the external 



or the internal perimeter, governs the behavior at Uc, and to decide whether the external 
perimeter dominates the interface for u > Uc, i.e. c^ep = di > dip, or if d-^p = di = dip in 
this regime. The values for c^ep and dip according to the above considerations are given in 
Table 



V. SATW IN ONE DIMENSION 

As we mentioned above, in = 1 there is no swelling- collapse transition, as the walk is 
collapsed for all u. The exponents k and u are k = 1/2 and u = 1/2 |TT| , |2T|] , in accordance 
with Eqs. (^. Following the considerations leading to Eq. (p!OD, the probability of the walker 
to be on the boundary of a S'y4rM^-cluster is {Ni) / {S) = 2/{S), therefore SATW in d = 1 
belongs to case (ii), since ai = Omax = oi = and = 1- Thus, = 1/(2 — Omax) = 1/2 as 



expected from analytical results |2T|. Note that, based on Eq. (H), one can derive (although 



not rigorously) a closed form expression for {S{t)) in ci = 1, extending the results of 



Since in ci = 1 the conditional probability to expand the cluster while being on a perimeter 
site is Pi = l/(exp(M) + 1) (cf. Eq. (0)) and A^i = 2, we obtain 

= _^ i_ (16) 

dt exp(u) + 1 {S{t)) ' ^ ^ 

yielding 

2t ^ 



This result is strongly supported by numerical simulations, see Fig. 
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In (i = 1 it is also possible to solve the inverse problem of deriving the average time 
{t{S)) to visit a fixed number of visited sites S, 

(((S)) = S-l+(^^i>Mii^. (18) 

2 Pi 

For asymptotically large S, Eq. ( pTsD yields {t{S)) ~ S'^, recovering the expected scaling 



with k = 1/2. To derive Eq. (|18|) we use the standard approach to RW in c? = 1 (see, for 
instance, Ref. ||3^): The number of sites S visited by a walker in d = 1 is equal to the span 
of the random walk. The time t{S) passed before the span reaches S can be represented as 
the sum 

s 

t{S) = J2x^, (19) 
i=2 

with Xi = t{i) — t[i — 1) being the time spent before the span increases from z — 1 to z, given 
that the walker is initially at the boundary. Thus, the mean time {t{S)) to visit S sites can 
be calculated as 

m)) = Y,^x^). (20) 

i=2 

with = and t(2) = 1. Let us consider the nature of the variable Xi some detail. 
When at the boundary, the walker increases the span with the probability of succeess Pi. 
With probability 1 — Pi the walker stays inside the cluster and undertakes an excursion 
until it hits the boundary again, which presents him with another opportunity to increase 
the span. We introduce the random variable /i which is the number of unsuccessful attempts 
before the span is increased. In other words, is the number of excursions into the cluster, 
and with probability Pi(l — Pi)™, m > 0, /x is equal to m. Now we can decompose the Xi ^ 



Xi 



where Tj is the lenght of the j-th excursion, and 1 is added to account for the jump from 
the boundary into the cluster before the excursion started. The random variable Tj can be 
viewed as the time a random walker started at site 1 with the boundaries at and i — 2 
needs to reach a boundary. Its mean is known to be (tj) = z — 3 Q]. The mean of Xi is just 
the average number of excursions multiplied by the average length of one excursion. More 
rigorously, we express (xi) in terms of conditional averages as follows 



H = m) P{fi = m) . (22) 



The variables tj are independent identically distributed random variables, therefore 



'J 

oo 



(x.) = l + ^m(r, + l)P(/i = m) = l + (^-2) ^ m A (1 - = 1 + ^^^^ 
1 1 Pi 

m=l m=l ^ 

(23) 

By combining Eqs. (^OD and (|^) one obtains Eq. ([TsD . 
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VI. CONCLUDING REMARKS 



In contrast to all other known random walk models with interaction, the SATW model 
exhibits a swelling-collapse transition at a critical attraction Mc in d > 2. The transition is 
similar to the swelling-collapse transition observed at the 'G point' T = O of SAW with 



an attraction term exp{—A/T), A < |jT2|,[T3|. It can only arise because the attractive 



interaction energy is of the same order as the configurational entropy in c? > 2. Below Uc the 
entropy dominates and the walk is in the universality class of simple RW, above Uc the energy 
governs the behavior and the walk collapses. At Uc, both contributions balance each other, 
leading to a new universality class. In d = 1, due to small number of possible configuration 
caused by the geometrical constraints, the walk is always collapsed, even without interaction 
u = 0. 

Analyzing the structure of the cluster grown by SATW in detail, we determined the 
fractal dimension of the cluster and its interface. In d = 2, the cluster is always compact 
while the interface has a fractal dimension di = 2 and 1.50±0.01 below and at u^, respectively, 
dominated by the internal perimeter. Above Uc, we found di = dip = c^ep = 1- The fractal 
dimension of the external perimeter at Uc, c^ep = 1-25 ± 0.05, is smaller than below Uc, 
dEP = 4/3, demonstrating once more the novel universality class at criticality. In d = 3 the 
cluster is compact for u < Uc, while for u > Uc the fractal dimension is df = 2. Probably the 
most unexpected result is that ind = 3 the interface has a fractal dimension di = 2 above and 
below Uc, whereas at criticality it increases to dj = 2.73±0.03. This appealing structure at Uc 
is of interest on its own regarding that many challenging problems in physics, chemistry, and 
biology are associated with growth patterns in clusters and solidification fronts. The results 
for the fractal dimensions of the perimeters could also be helpful when investigating the 
structure of the cluster grown by ^^VF at the corresponding transition, where sufficiently 
large clusters are difficult to simulate because of the attrition of the walkers when stepping 
into its own dead ends. 
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TABLES 





U < Uc 


U = Uc 


U > Uc 


V 

d = 2 k 


1/2 
1 


0.40 ±0.01 
0.80 ±0.01 
Uc = 0.88 ± 0.0. 


1/3 
2/3 


u 

d = 3 k 


1/2 
1 


0.303 ± 0.005 
0.91 ±0.01 

Uc = 1.92 ± o.o; 


1/4 
3/4 



TABLE I. The exponents v and k as well as the estimated values for the transition point Uc 
for SATW in d = 2 and d = 3. 









U < Uc 


U = Uc 


U> Uc 
















at 




d = 


2 


l<i <2d 
i = 2d 


0.95 (1) 
1 


0.90 (1) 
0.95 (1) 


0.75 
1 


0.60 
0.80 


0.50 
1 


0.33 
0.67 


d = 


3 


l<i <2d 
i = 2d 


1 
1 


1 
1 


0.91 

1 


0.83 
0.91 


0.67 
1 


0.50 
0.75 



TABLE IL The exponents Oj and k^ for 1 < i < 2d and i = 2d for all three regimes u < Uc, 
u = Uc and u > Uc in d = 2 and d = 3 from numerical simulations. The errors are of the order 
of one percent. For u < Uc in d = 2, most probably due to logarithmic corrections it is not 
possible to obtain the asymptotic value of the exponents in numerical simulations, therefore we give 
them in parentheses. Note that the results for a for the regimes u < Uc and u > Uc are supported 
by combining k = 1/2 and Eq. ( |To| ) as well as Eqs. (|^) and ([lO|), respectively. 





U < Uc 


U = Uc 


U > Uc 




df 


2 


2 


2 


d = 2 


di 


2 


1.50 ±0.01 


1 


cZep 


4/3 


1.25 ±0.05 


1 




dip 


2 


1.50 ±0.01 


1 


d = 3 


di 


2 


3 


3 








di 


2 


2.73 ±0.03 


2 



TABLE III. The fractal dimension d{ of the cluster and di = ai df of the interface grown by 
S'ATW for all three regimes u < Uc, u = Uc, and u > Uc in d = 2 and d = 3. Additionally 
in d = 2 the fractal dimensions of the external perimeter c^ep and the internal perimeter dip are 
given. In d = 3 it is particularly not clear whether the external or the internal perimeter dominates 
the interface at u = Uc- In case no errorbars are given, the numerical results are supported by 
analytical considerations. 
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FIGURES 



FIG. 1. Representative example of the structures grown by SATW after t = 2 ,2 , and 2 
time steps for u = 0, u = 0.88 (= Uc), and u = 2.5 in d = 2. Gray sites form the interface, black 
sites are completely surrounded by already visited neighbor sites. The length scale is arbitrary 
chosen so that the clusters fill the size of the box. Note that for u < Uc, one can easily follow the 
growth process because of the distinctive structure of the cluster, while for u = Uc it is difficult to 
follow, since the walker keeps coming back more often altering the structure. For u > Uc it is not 
possible to do so as the grown clusters are compact. 



FIG. 2. The values of the exponents k (squares) and (circles) vs attraction u in d = 2 (open 
symbols) and d = 3 (filled symbols) for t = 10^, obtained by a least square fit of the slope of 
ln(i?^(t)) and ln(S'(i)) vs Int for large t, respectively. Note that for u > Uc the values of k and v 
approach the theoretical predictions of Eqs. (§), marked as dashed lines. The extimated values of 
Uc are = 0.88 ± 0.05 in d = 2 and = 1.92 ± 0.03 in d = 3. 



FIG. 3. The average number of cluster sites {Ni{t)) having i already visited nearest neighbor 
sites plotted vs the average number of all cluster sites {S(t)) up to t = 2 • 10^ time steps (for 
u < Uc in d = 2 and d = 3 only up to t = 10^ and t = 3 ■ lO"^, respectively) averaged over 100 
configurations in d = 2 for (a) u = 0.5 < Uc, (b) u = 0.88 = n^, (c) u = 2.5 > Uc, and in d = 3 for 
(d) u = 1 < Uc, (e) n = 1.92 = Uc, and (f) n = 4 > Uc. In d = 2 the data for the different values 
of i are marked by i = 1 (circles), i = 2 (squares), i = 3 (diamonds), and i = 4 (pluses), whereas in 
d = 3, i = 1 (circles), i = 2 (squares), i = 3 (diamonds), 1 = 4 (upward triangles), i = 5 (downward 
triangles), and i = 6 (pluses). Note that (d) clearly differs from (a),(b),(c), and (f). The value of 
the slopes determined from the data are summarised in Table 



FIG. 4. (a) The successive slopes = dln(A'^j)/dln(S') of the data from Fig. |^(e) plotted 
vs l/ln(5'). A linear extrapolation of the points to the limit l/ln(S') yields our estimates 
Oj = 0.91 lb 0.01 for 1 < i < 2d and a2d = 1-00 it 0.01, clearly revealing that this regime belongs 
to case (ii). The data is marked as in Fig. ^(e). (b) The average number of cluster sites having i 
already visited nearest neighbors divided by the average number of all cluster sites {Ni{t)) / {S{t)) 
vs {S{t)) of the data from Fig. |3|(a). Note that the data for i = 4 marked by pluses shows a 
dominant behavior compared to the data for z = 1 (circles), i = 2 (squares), and i = 3 (diamonds), 
implying that the regime u < Uc in d = 2 belongs to case (ii). 



FIG. 5. The successive slopes dj = dln{I{R))/d\nR of the mean mass {I{R)) of the interface 
inside a disk resp. a sphere of radius R plotted vs 1/R for u = Uc at t = 3 - 10^ resp. t = 3 - 10^ time 
steps averaged over 1000 configurations. A linear extrapolation of the points to the limit 1/R 
yields our estimates (a) di = 1.50 it 0.02 in d = 2 and (b) di = 2.72 it 0.04 in d = 3, confirming 
the results summarized in Table III, which are marked by an arrow. 



14 



FIG. 6. The successive slopes oep = dln(A^Ep)/dln(5') plotted vs l/ln(5) averaged over 500 
configurations for u = 0.5 < Uc (diamonds) up to t = 10^ time steps, for u = Uc (squares) up to 
t = 10^ time steps, and u = 2.5 > Uc (circles) up to t = 2 • 10^ time steps in d = 2 . A linear 
extrapolation of the points to the limit l/ln(S') yields our estimates oep = 0.68 ± 0.02 and 
0.50 lb 0.02 for below and above Uc, respectively, confirming the results summarized in Table ID, 
which are marked by an arrow. At Uq, oep is clearly smaller than for below Uc- 



FIG. 7. The successive slopes for c^ep = dln(A^Ep)/dlni? and dip = dln(A^ip)/dlni? of the 
mean mass (A^ep) and (Ajp) of the external and internal perimeter, respectively, plotted vs 
for (a) u = Uc at t = 10* time steps and (b) u = 2.5 > Uc sd, t = 2 ■ 10^ time steps averaged 
over 500 configurations in d = 2. A linear extrapolation of the points to the limit 1/R 
yields our estimates (a) dip = 1.48 ± 0.05 (circles) and (b) c^ep = 1.02 it 0.03 (squares) and 
dip = 1.01 lb 0.03 (circles), confirming the results summarized in Table ID, which are marked by 
an arrow. At Uc, (Iep (squares) is determined to be dEP = 1-25 ± 0.05. 



FIG. 8. The average number of visited sites {S{t)) in d = 1 for different values u, u = 2 (cir- 
cles), M = 4 (squares), u = 6 (diamonds), and u = 9 (triangles), plotted vs time t scaled by 
(1 + exp(u))^ shows a good data collapse in agreement with Eq. ( p^ (solid line). The plateau for 
small t corresponds to the average time (equal to 1 + exp(u)) needed for the walker to escape the 
initial cluster of the size 2. 
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